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ABSTRACT 


<10 

This  paper  provides  a basic  description  of  the  operation 
of  a shearing  interferometer  used  for  phasefront  sensing.  The 
signal  processing  of  the  interferometer's  outputs  are  studied. 

A basic  demodulation  scheme  and  a phaselocked  loop  configuration 
are  discussed  and  compared.  Some  elementary  phaselocked  loop 
theory  is  provided.  A possible  configuration  using  optical 
feedback  for  determining  off-axis  tilt  of  a plane  is  presented. 
The  signal- to-noise  ratio  of  the  shearing  interferometer's 
outputs  is  improved  by  a slight  modif ication  of  the  optical 
configuration.  Finally,  a geometric  representation  of  a 
distorted  phase  profile  is  provided  which  will  describe  the 
deleterious  effects  of  a planar  phasefront  propagating  through 
the  turbulent  atmosphere. 
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I Introduction 

The  turbulent  atmosphere,  which  's  sometimes  referred  to 
as  a random  inhomogeneous  medium,  distorts  and  deforms  optical 
radiation  as  it  propagates  outward  from  its  source  (Refs.  1-3). 
That  is,  if  a source  is  radiating  quasi-monochromatic  light 
(Ref.  4,  pg  13),  a distorted  plane  wave  will  be  received  in  the 
far  field.  If  this  source  were  a laser,  then  its  far  field 
target  would  not  receive  a focused,  intense  beam  because  of 
atmospheric  turbulence.  The  problem  of  maximizing  the  intensity 
of  a laser  beam  on  a far  field  target  is  evident. 

It  has  been  shown  through  reciprocity  (Ref.  5)  that  an 
optical  system  which  projects  a laser  beam  onto  a far  field  target 
will  experience  identical  atmospheric  turbulence  effects  as 
the  image  of  a radiating  target  on  the  optical  system.  There- 
fore, if  one  can  measure  the  phase  profile  of  a target's 
radiation,  then  the  effects  of  the  turbulent  atmosphere  will  be 
known.  Compensation  for  these  effects  can  then  be  made  so  that 
a laser  beam  will  not  suffer  the  same  deleterious  effects  but 
will  instead  have  maximum  intensity  on  the  target. 

A common  configuration  for  the  optical  system  which 
compensates  for  atmospheric  turbulence  is  shown  in  Figure  1. 

The  telescope  acts  as  the  exit  aperature  for  the  laser  beam 
and  the  receiving  aperature  of  the  target  radiation.  The  optical 
radiation  from  the  target  enters  the  telescope,  reflects  off  the 
tilt  mirror,  reflects  off  the  deformable  mirror,  and  upon 
incidence  to  the  beamsplitter  is  reflected  into  the  phasefront 
sensor.  The  phasefront  sensor  reconstructs  the  phase  profile 
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of  the  target's  radiation.  Through  appropriate  signal  processing 
of  the  reconstruction,  the  atmospheric  turbulence  can  be 
compensated  for  so  that  the  intensity  of  the  laser  beam  is 
maximized  on  the  target.  This  compensation  is  accomplished  by 
control  signals  to  the  tilt,  focus,  and  deformable  mirrors. 

The  laser  beam  propagates  through  the  beamsplitter,  reflects 
off  the  deformable  mirror,  reflects  off  the  tilt  mirror,  and 
exits  the  telescope  to  propagate  through  the  turbulent 
atmosphere  to  the  target. 

The  device  of  concern  in  this  paper  is  the  phasefront 
sensor.  Specifically,  a type  of  phasefront  sensor  referred  to 
as  a shearing  interferometer  is  studied.  Also  the  signal 
processing  of  the  interf eromeier  output  is  analyzed. 

The  following  chapter  provides  a description  of  the 
shearing  interferometer.  Chapter  3 discusses  the  processing 
of  detector  signals  using  Phaselocked  loop  and  Optimal  Linear 
Filter  theory.  A geometric  representation  of  a distorted 
phase  profile  is  discussed  in  Chapter  4. 
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II  Phasefront  Measuring  by  use  of  a Shearing  Interferometer 

There  are  no  direct  detectors  that  can  measure  the  phase 
of  a field.  Detectors  measure  the  intensity  of  a field  and  thus 
phase  information  is  .lost.  The  shearing  interferometer  provides 
a solution  to  this  problem. 

A shearing  interferometer  (Refs.  6-9)  is  a phasefront 
sensor  whose  outputs  can  be  used  to  reconstruct  the  phase 
profile  or  contour  of  the  input  optical  field. 

The  basic  components  of  a shearing  interferometer  are  two 
Fourier  transform  lenses,  a rotating  Ronchi  grating  (Ref.  10, 
pg  207),  and  a detector  array  as  shown  in  Figure  2. 

In  Figure  2 Zq,  z^,  and  are  specific  locations  along 
the  z or  optic  axis,  and  f is  the  focal  length  of  lenses 
and  L^-  The  Ronchi  grating  consists  of  a circular  disk  with 
alternating  clear  and  opaque  radial  lines  extending  from  near 
the  center  to  the  outer  edge  (Ref.  9)  as  shown  in  Figure  3, 

pg  6. 

The  rotating  Ronchi  grating  which  has  radius  r,  rotating 
at  Wq  radians  per  second,  and  having  N pairs  of  radial  spokes 
will  have  a linear  velocity  of 

v = wQr  (1) 

and  a grating  period  of 

d = 2Trr/N  (2) 

At  any  radius  the  amplitude  transmittance  function  of  the 
grating  can  be  approximated  by  a square  wave,  which  can  be 

If 




r 

R 

g 


represented  as  the  following  Fourier  series 

g(x^-vt)  = E An  exp  [+j?irn(x^-vt)/d]  (3) 

where  the  A 's  are  Fourier  coefficients.  The  coefficients 
n 

for  the  Ronchi  grating  are  shown  in  Table  1.  (Ref.  11,  pg  28) 
In  the  following  analysis  constant  phase  delay  terms  have 

5 


i 


Figure  3-  Rotating  Ronchi  grating 

been  neglected.  A time  dependence  of  exp[j27rfQt]  for  the 
incident  quasimonochromatic  field  is  assumed.  When  an  incoming 
field  w(x0,yQ)  = A±nc  ( xQ  ,y0  )exp  [+ j <}>  (xQ  5yQ  ) ] (where  Ainc  is  a 
real  amplitude  and  XqjYq  are  coordinates  at  Zq  of  the  optic 
axis)  is  located  at  f in  front  of  lens  L^,  the  Fourier 


order  n 

amplitude  of  An 

phase  of  A^ 

intensity 

0 

.5000 

0° 

.2500 

±1 

1/tt 

±90° 

.1013 

±2 

0 

- 

- 

±3 

1/3* 

±90° 

.0112 

r 


Table  1.  Fourier  Coefficients  for  a Ronchi  gratin, 


transform  of  the  wavefront 


F [w ( x0 , yQ  ) ] 


W(x1,y1) 

x1/Xf 

yx/xf 


(4) 


is  proportional  to  the  field  at  the  plane  of  the  rotating 
grating.  (x-j^  and  are  coordinates  at  of  the  optic  axis, 
fx  and  f are  spatial  frequencies  , and  X is  the  optical 
wavelength)  Immediately  following  the  grating  the  field  becomes 


UU-j^Yj)  a WU^y^-Z  Anexp  [ j2irn(x1-vt)/d]  (5) 

At  the  detector*  plane  (at  location  z2  ° f °ptic  axis)  the 
field  becomes 

F[U(x1,y1)]  * w(-x2,-y2) 

fx  ■ X2/Xf 

fy  - y2/»f 

© F[g(x1-vt)]  (6) 

fx  = x2/Xf 

fy  - y2/xf 

where  the  symbol  © represents  a two  dimensional  convolution 
(Ref.  11,  pg  82).  The  Fourier  transform  of  equation  (3)  becomes 


F [gtx^-vt ) ] “ G(x2-vt) 

fx  ■ X2/ir 

fy  - y2/u 

= I exp[-j(k/f )x2vt]An 

• 6(x2-nAf/d)S(f  ) (7) 

«y 

The  above  equation  is  nonzero  only  when  x2  = nAf/d,  and  by 
substituting  for  x2  in  the  first  term  of  equation  (7)  we  have 
exp  [-j  (k/f  )x2vt  ] = exp  [-j  2irnvt/d]  where  wQ  = 2imv/d  is  the 
modulated  frequency  of  the  electric  field.  Now  equation  (6) 
becomes 

U(x2,y2)  = FlUCx-^y^]  w(-x2,-y2) 

fx  ■ X2/Af 

fy  - y2/xf 

@ l exp[-jwnt] 
n=0 ,+1,-1 

• A 6 (x0-nAf/d) 

n 2 

• 6(f)  (8) 

«/ 


where  n is  restricted  to  0,  +1,  and  -1  since  higher  orders 
result  in  higher  frequencies  which  are  filtered  out  by  the 


where  X is  the  wavelength  of  the  field.  Figure  4 shows  the 


shearing  pattern  of  the  first  three  terms  produced  at  the  detector 
plane.  The  three  electric  fields  of  Figure  4 are 

v0  = A0Ainc(x2Jy2)exp[j<t>(x2,y2)] 

V+1  = A+1Ainc(x2  + s,y2)exp[j<J)(x2+s,y2)  ]exp[-jw3t] 

and 

V_1  = A_1Ainc(x2-s,y2)exptj4.(x2-sJy2)  ]exp[+jwQt]  (10) 

where  Aq,  A+1,  and  A ^ are  found  in  Table  1.  The  total  electric 
field  at  the  detector  plane  is 

Ud^x2,y2^  = V0  + V+1  + V-1  + ^hiSher  or,der  terms}  (11) 

where  the  higher  order  terms  are  neglected  as  they  will  be 
filtered  out  by  the  electronics.  Since  the  detectors  of  the 
shearing  interferometer  measure  intensity,  the  intensity  of 
Ud(x2,y2)  becomes 
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■'"d^X2,y2^  ~ ^x2  * y2  ; ^d  ^ x2  ’ y2  ^ 


<V0  + V+1  + v-i][vo  + V+1  + V-l^ 


y y*  + V v*  + V V*  + V V*  + V V* 
v0  0 v+lv+l  -1-1  v+lv0  +1  -1 


+ V.iVg  + + v0v*1  + VqV*,  (12) 


where  * denotes  the  complex  conjugate  of  the  designated  quantity. 

For  simplicity  let  Ainc(-X2,y2^  = Alnc (*2+s ,y2 ) = Ainc (x2~s ,y2 ) 

= A^,  then  equation  (12)  becomes 


2A 

Id^x2,y2^  = At  + ~T~ 


+ — cos(wnt 

IT  U 


+ .time  varying  terms  of, 
* frequency  2Wq 


- 90°  + <j>(x2,y2)  - 4> (x2+s , y2 ) ) 


A 


+ cos  ( Wgt  - 90°  + 4>  ( x2-s , y2 ) - <t>  ( x2 , y2 ) ) 

(13) 


(See  Appendix  for  result  when  the  sheared  amplitudes  are  not 

equal.)  The  DC  term  A^  + 2A|/tt  , and  the  double  frequency  terms 

2wgt  are  also  filtered  out  by  the  electronics,  thus  equation  (13) 

reduces  to  p , . , . 

2 A.  4>  (x?-s,yP)  - (Xp+s  ,yn  ) 

Ia(x2-y2)  ■ - 90‘  + — 5 -> 


<j)(x2-s,y2)  + <j>(x2+s,y2) 
c o s ( — 

(14) 


Assuming  the  phase  variation  across  a shear  distance  is  small 
(that  is  U(x2-s,y2)  + ( x2  + s ,y?  ) )/2  = 4>  ( x2 , y2 ) ) equation 
(14)  becomes 
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2 A.  <}> (x„  — s , y P ) - |(x.+s,y5) 

Id(x2,y2)  = 7^'sin(wOt  + 2 (15) 

This  result  shows  that  the  overlapping  fields  at  the  detector  plane 

produces  a signal  whose  phase  is  proportional  to  the  phase 

difference  or  slope  of  the  wavefront.  Therefore,  the  shearing 

interferometer  measures  a set  of  phasefront  profiles  in  the 

direction  of  the  shear. 

A phasefront  sensor  uses  two  shearing  interferometers  that 
produce  shear  in  orthogonal  directions  (x  and  y)  providing  a 
two-dimensional  representation  of  a distorted  phasefront.  For 
simplicity  assume  that  there  are  twelve  detectors  (six  for  the  x- 
channel  interferometer  and  six  for  the  y-channel  interferometer) 
so  that  the  phasefront  sensor  outputs  twelve  signals  whose 
phases  are  phase  differences  that  define  the  slope  of  the  phase- 
front  in  two  dimensions  (Ref.  6,  pg  668).  In  Figure  5 the  detector 
arrangements  are  shown  along  with  a phase  map  illustrating  how  the 
phase  difference  at  each  detector  is  derived.  Detectors  are 
located  behind  each  subaperture  labeled  for  the  x and  y 

shear. 

In  order  to  obtain  the  actual  phasefront  a reconstruction 
process  is  necessary  which  uses  the  twelve  phase  difference 
measurements  as  shown  by  Hardy  (Ref.  6,  pg  669).  The  reconstructed 
phasefront  can  be  represented  by  a set  of  orthonormal  polynomials 
which  will  be  described  in  Chapter  4.  These  polynomials  will 
describe  the  amount  of  random  tilt  (x  and  y),  defocus,  and 
astigmatism  that  the  received  phasefront  has  suffered  through 
its  propagation  through  the  turbulent  atmosphere. 
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Figure  5.  ?hase  difference  measurements  and  location  of  detectors 


Since  it  is  the  phase  that  is  of  interest,  equation  (15) 
for  each  of  the  twelve  detectors  must  be  demodulated  in  such 
a way  that  the  phase  term 


* 
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(16) 


<t>(x2-s,y2)  - <j>(x2+s,y2) 

2 

is  determined. 

The  actual  intensity  at  the  detector  plane  is  very  small 
compared  to  that  at  the  input  of  the  shearing  interf erometer 
due  to  the  less  than  unity  transmittance  of  the  optics. 

It  can  be  shown  that  employing  a sinusoidal  grating 
instead  of  the  Ronchi  grating  can  increase  the  intensity  at  the 
detector  plane.  By  using  appropriate  Fourier  coefficients  and 
equation  (.3),  a result  for  any  grating  can  be  obtained. 

The  next  chapter  discusses  possible  configurations  for 
determining  the  phase  of  a detector  signal. 


Ill  Processing  of  Detector  Signals 


Basic  Demodulation  Scheme 

A commonly  proposed  system  for  determing  the  phase 
of  a signal  is  shown  in  block  diagram  form  in  Figure  6.  The 
detector  signal  s(t)  = DsinCwgt  + e(t))  (where  D is  a real 
amplitude  and  6 ( t ) = [ $ ( t , x^-s  ^ ) ~ <t>  ( t , x^+s  ,y^ ) ] /2  ) 

enters  a bandpass  filter  of  center  frequency  Wq/2tt.  Then 
the  signal  is  synchronously  demodulated  into  two  quadrature 


Figure  6.  Basic  demodulation  scheme  for  detector  signals 


components  Dsine(t)  and  Dcos9(t),  and  lowpass  filtered  to 
eliminate  some  of  the  noise.  The  A/D  ROM  (analog-to-digital 
read-only-memory)  performs  an  inverse  tangent  operation  through 
a "table  look-up"  to  determine  0. 

For  simplicity  of  analysis  this  scheme  is  studied  on  a 
continuum  in  time.  Assuming  the  presence  of  bandpass,  zero 
mean  noise,  the  actual  detector  output  becomes 


s(t)  = Psin(wQt  + 0 ( t ) ) + n1(t)sinwQt 

+ n2(t)coswQt  (17) 

where  n^(t)  and  n2(t)  are  statistically  independent  zero  mean 
random  processes.  The  output  of  the  scheme  of  Figure  6 would 
be  (Ref.  12,  pg  174-178) 

y(t)  = e(t)  + tan~1(^-jt^,^tj)  (18) 

where 

n'(t)  = n^(t)cose(t)  + n0(t)sine(t) 

n"(t)  = -n1(t )sin0(t ) + n2(t)cos0(t)  (19) 

and  n'(t)  and  n"(t)  are  zero  mean  random  processes  with  distri- 
butions equal  to  those  of  n^(t)  and  n 2 ( t ) . For  the  high  signal 
to  noise  case  (that  is,  Prob{n'>D}  = Prob{n">D}  - 0)  then 
equation  (18)  can  be  approximated  by 

y(t)  = 0 ( t ) + (20) 

The  variance  of  equation  (20)  is  derived  as  follows 

E [ y ( t ) 1 = e(t)  (21) 


E[y2(t)]  = E[02(t)  + 2 6 ( t ) 


-n, (t ) sin0 (t ) + n-(t)cos0(t) 


(- 


D 


-n, (t )sinfl (t ) + no(t)cos0(t)  2 

+ <— D > 1 


15 


where 


= 02(t)  + i~E[n2(t) 
D 1 

= 02(t)  + o2/D2 


E [n  ] = Rn(0) 


-Bl/2  N B 

= f S (f)df  = -3—^ 
-B.,/2  ^ 


and  B-^  is  the  bandwidth  determined  by  the  lowpass  filter  of 
Figure  6.  Therefore,  the  variance  of  y(t)  becomes 


2 NnB, 
0 = 0 1 

y — 2~ 

2D 


This  variance  is  an  indication  of  how  much  on  the  average  the 
random  process  y(t)  fluctuates  about  its  mean  0 ( t ) . The  band- 
width B-^  must  be  selected  large  as  compared  to  the  bandwidth  of 
a scheme  which  uses  a pX-ioKl  knowledge  of  the  process  concerned. 

The  phaselocked  loop  (PLL)  is  a scheme  that  uses  some  a pfiloxl 
information  and  is  discussed  in  the  next  section.  Some  basic 
PLL  theory  will  be  provided  first  (Refs.  13,  pg  66  — 8-4 ; 1*J,  pg  127-131) 

Basic  Phaselocked  Loop  Theory 

The  Phaselocked  loop  configuration  is  used  for  measuring 
the  phase  of  a signal.  The  basic  elements  of  a PLL  are  a 
phase  detector,  a filter,  and  a voltage  controlled  oscillator 
(VCO)  as  shown  in  Figure  7.  (The  phase  detector  is  represented 
as  a multiplier . ) 


The  effective  operation  is  dependent  on  the  amount  of  a pK.lon.1 
information  that  can  be  designed  into  the  circuit.  The  PLL 
provides  an  estimate,  r(t)  of  the  input  signal  s(t),  and  by 
comparison,  the  loop  continously  modifies  the  local  reference 
(VCO  output)  to  reflect  changes  in  the  input.  An  error  signal 
e(t)  is  produced  at  the  output  of  the  phase  detector  (the 
multiplier  of  Figure  7)  if  the  input  signal  leads  or  lags  the 
local  reference.  The  error  signal  is  filtered  and  sent  back  to 
the  VCO  in  such  a way  as  to  change  the  frequency  of  the  VCO  so 
that  it  approaches  the  frequency  of  the  input  signal. 


r 


In  order  to  determine  the  loop  equation  of  the  system, 
s(t)  is  considered  to  be  a narrowband  random  process 


s(t  ) 

= Dsin(wQt  + 6 (t ) ) 

(25) 

r(t) 

= K^cos(wQt  + 0 ( t ) ) 

(26) 

where  6(t)  is  the  loop  estimate  of  the  input  signal  phase  e(t) 
in  the  absence  of  noise,  and  is  the  amplitude  of  the  VC 0 

output.  The  output  of  the  multiplier  e(t)  neglecting  double 
frequency  terms  becomes 

e = DK1Msin(0-0)  (27) 

where  M is  the  multiplier  gain.  The  argument  (6-6)  is  called 
the  total  loop  phase  error.  The  output  of  the  filter  with  time 
response  f(t)  becomes 

z(t)  = e(t)F(s)  (28) 

where  F(s)  is  the  Laplace  transform  of  f(t)  and  the  notation 
used  is  consistent  with  that  of  Lindsey  (Ref.  13).  (That  is, 
equation  (28)  is  equivalent  to  z(t)  = e(t)  © f (t ) . ) The  output 
phase  of  the  VCO  6 referenced  to  zero  is  regarded  to  be  a linear 
function  of  the  control  signal  z(t).  When  z(t)  is  zero  (which 
implies  that  e(t)  is  also  zero)  then  the  VCO  outputs  a constant 
quiescent  frequency  of  Wq  radians/second.  As  mentioned  earlier 
Wq  is  also  the  frequency  of  the  signal  output  from  a detector. 
When  z(t)  is  nonzero,  then  the  VCO  output  becomes  Wq  + K^^(-)z(t) 
where  is  the  gain  of  the  voltage  controlled  oscillator 

(VCO)  in  units  of  radians  per  second  per  volt.  The  phase 
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estimate  9 developed  by  the  loop  is 


0 


KVCOZ 


s 


(29) 


where  again  the  notation  used  is  consistent  with  that  of 
Lindsey.  (The  1/s  term  is  Laplace  notation  for  the  operation 
of  integration.)  Substituting  for  equation  (27)  and  (28)  into 
(29)  yields 


9 = 


DK1KVC0MP(s) 


sin( 9-9 ) 


= 5Milsin$  (30) 

s 

where  K = is  the  open  loop  gain  and  $ = (0-0).  Since 

9 = 0-$,  equation  (30)  becomes 


$ 


_ DKP(s)  . . 

0 - — -sin$ 

s 


(31) 


Equation  (31)  is  the  equation  for  the  PLL  in  the  absence  of 
noise . 

In  the  presence  of  noise,  the  input  of  Figure  7 becomes 
x(t)  = s(t)  + n^(t).  Now  the  error  signal  e(t)  depends  not 
only  on  the  total  loop  error  <f(t)  but  also  the  noise.  The  noise 
n^(t)  is  assumed  to  be  a narrowband  random  process 

n^(t)  = nccosWgt  - ngsinwgt  (32) 

= Ni(t)cosCwQt  + 0^(t))  (33) 


where  N.(t)  is  the  envelope  of  the  noise  and  0<(t)  is  the  random 
1 


19 


\ 

A 


phase.  From  equation  (33)  we  have 

n = N,cos0. 
c i 1 

and 

n = N.cosO. 
s i i 

The  input  now  becomes 


(34) 


x(t)  = Dsin(wQt  + e(t) ) + N1(t)cos(wQt  + Q^(t)) 

= (D-Ns ) sin(wgt  + 0(t))  + Nccos(wQt  + 0(t))  (35) 

where  Nc  = N^cos(@i-0)  and  Ns  = N^sin(0^-0).  The  error  signal 
then  becomes 

e = K, M [ (D-N  ) sin( 0-0 ) + N cos(0-0)]  (36) 

_L  S C 

Once  again  the  double  frequency  terms  have  been  neglected. 

Since  0(t)  = K^^z  ( t )/s , and  z(t)  = F(s)e(t),  then 

9 = KF^s) [Dsin( 0-0 ) + N ( t ) ] (37) 

where 

N(t)  = N cos(0-0)  - N sin(0-0)  (38) 

C 3 

From  equation  (38)  it  is  obvious  that  the  noise  is  a function 
of  the  phase  error  (0-0).  Rearranging  equation  (37)  the 
following  is  a nonlinear  integrodif ferential  equation  of  the 
system 

(0-0)  = $ = 0 - --F(.s  )-[Dsin(  0-0  ) + N ( t ) ] (39) 

s 

The  baseband  loop  model  which  equation  (39)  represents  is 
shown  in  Figure  8. 
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If  the  phaselocked  loop  is  capable  of  reducing  the  phase 
error  to  a small  value,  then  the  approximation  sin$  - 4>  becomes 
valid.  Equation  (39)  then  becomes 

$ & (0-e)  = G - s ^ [D$  + N cos$  - N sin#] 

s c s 

0 _ KF(s)  [D*  + N ] (40) 

S 

since  sin$  = $ implies  that  cost  = 1,  and  N >>  N $ . The  linear 
baseband  model  which  represents  equation  ( 4 0 ) is  shown  in 
Figure  9- 


Figure  9.  Linear  baseband  PLL  model  in  the  presence  of  noise 

It  can  be  shown  (Ref.  13,  pg  81-82)  that 

N 

Rn  (t)  ■ 6(0  (41) 

c 

Therefore,  the  noise  process  { Nc(t,$)  } is  white.  Solving 
equation  (40)  for  $ yields 

^ (42) 

D 

The  closed  loop  transfer  function  of  Figure  9 is 
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„/  >.  A0(s)_  DKF(s) 

¥tsT  s + DKP(s) 


Therefore  equation  (42)  becomes 


(43) 


4>  = closed  loop  phase  error 

Nc 

- I1  - Ve  ~ 


(44) 


where  phase  error  can  be  caused  by  both  the  process  0 
changing  faster  than  the  loop  can  track  and  by  the  noise  N . 
Equation  (44)  shows  that  there  will  be  no  contribution  to  the 
phase  error  by  0 as  long  as  [1  - H$ ] is  sufficiently  highpass 
so  that  its  spectrum  doesn't  overlap  with  that  of  0.  For  this 
case  the  phase  error  will  be  dependent  solely  on  the  noise. 

The  variance  of  the  phase  error  in  this  case  is 


+B„/2 

/ 

-B2/2 


iVs> 


, SN  (s) 
2 c 

D2 


ds 


BpNn 

= = 1/SNR  (45) 

2D 

where  B^  is  the  loop  bandwidth,  and  N^/2  is  the  power  spectral 

density  of  N . As  seen  in  equation  (45)  the  phase  error  (due 
c 

to  noise  only)  is  inversely  proportional  to  the  signal  to 
noise  ratio.  The  loop  bandwidth  B2  should  be  considerably  less 
than  of  equation  (24)  since  some  a pHA.oH.-L  statistics  have 
been  designed  into  the  PLL  circuit. 

The  phaselocked  loop  has  been  shown  to  exhibit  the  property 
of  threshold  extension,  that  is  smaller  signals  can  result  in 
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similar  performance  (Ref.  15,  pg  359)-  This  extension  is 
typically  on  the  order  of  2 or  3dB  (Ref.  14,  pg  291). 

In  the  next  section  the  filter  F(s)  of  Figure  9 is  designed 
so  that  the  resulting  variance  of  the  phase  error  is  minimized. 
The  section  will  include  the  contributions  of  both  6 and  noise  to 
the  phase  error. 


Optimal  Filter  Design  for  the  PLL  Scheme 

In  this  section  the  filter  F(s)  of  Figure  9 is  designed 
so  that  the  resulting  variance  of  the  phase  error  is  minimized. 
If  $ of  equation  (44)  is  a random  process  then  (letting  E[0]=O) 

E[*]  = 0 (46) 


since  is  a zero  mean  white  process  with  autocorrelation 


N 


Rn  (t)  = E[N  (t)N  (t  + T)]  = -^S(t) 

C 


(47) 


From  equation  (44)  we  have 


oo 

E[*2]  = o2  = |1  - H$(s)|2SQ(s)ds 


+ 00 


+ f_m  |h4(s)|  Sn,(s)ds 


(48) 


where 


Sn'(s)  = 


N.  SN  <s) 

0 C 

2 ” 2 
2D  D 


(49) 


Equation  (48)  is  to  be  minimized  by  the  following  procedure. 

The  random  process  0 is  assumed  to  be  zero  mean  with  a 
first  order  Butterworth  power  spectral  density  (Ref.  16,  pg  502) 
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r 


2Q/ar 


s0 (s ) = 

6 1 + (s/a0)‘ 


(50) 


where  Rq(t)  = F[Sn(s)]  = Qexp[-aQT],  aQ  is  the  bandwidth  and 


'6"'  ‘ 1 e 

the  variance  of  d is 


°9  - V0)  ’ « 


(51) 


With  the  above  statistics,  a simple  filter  model  for  the  random 
process  0 is  derived.  Let  w(t)  be  white  gaussian  noise  with 


and 


E [w(t ) ] = 0 , 


E [w(t )w(t+T  ) ] - Rw(x)  = 2a0Q6(t) 


Sw(s)  = 2a0Q 


(52) 

(53) 
(5*0 


There  exists  a filter  H(s)  such  that  the  following  is  true 


and 


SQ(s)  = |H(s)rSw(s) 


H(s)H«(s)  = S(s)/S(s) 
0 w 


(55) 

(56) 


where  * denotes  the  complex  conjugate  of  the  designated 
quantity.  Substituting  into  equation  (56)  gives 

2 


H( s )H* ( s ) = 


1/a 


0 


1 + (s/aQ) 


2 


(57) 


Therefore 


H(s)  = 


s + a. 


(58) 
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Figure  10.  Model  of  the  random  process  8 

The  process  0 is  therefore  modeled  as  shown  in  Figure  10. 

From  Figure  10  we  have 

e(s)  = l (59) 

wTsT  s + aQ 

and 

e(s)  = -aQe(s)  + w(s)  (60) 

which  is  the  first  order  state  variable  equation  defining  the 
random  process  0. 

2 

In  order  to  find  the  minimum  value  of  a,  , continuous 

$ ’ 

Kalman  filter  equations  are  used  to  design  an  optimal  filter 
F(s)  in  Figure  9.  (See  Ref.  17,  pg  103-12^  for  the  development 
of  the  continuous  Kalman  filter  equations.)  From  Figure  9 we 
have 


L 


2 6 


(61) 


y(t)  = D0(t)  - D0(t)  + Nc(t,®) 

Using  Table  4.3-1,  pg  123,  Ref.  17,  the  state  estimate  becomes 

© = -a09(t)  + K2(t)[y(t)  - ( D 6 ( t ) - D0(t))] 

= -aQ0(t)  + K2(t)y(t)  (62) 

where  K^(t)  = 2P(t)D/NQ  is  the  Kalman  gain  and  P(t)  is  the 
variance  of  $ . The  Ricatti  equation  for  P(t)  becomes 

P(t)  = -2aQP(t)  + 2a0Q  - -2P-^it  )D-2-  (63) 

For  steady  state  set  P(t)  = 0 and  the  minimum  value  of  the 
variance  becomes 


a2  = P , . 

$ steady 

state 


N0a0 

2D2 


2D 


N0a0 (N0a0 


4D2Q) 


(64) 


where  the  plus  sign  of  the  quadratic 
since  the  variance  must  be  positive, 
gain  becomes 


2steady 

state 


2P 


steady^* 

state 


equation  has  been  chosen 
The  steady  state  Kalman 


-a. 


Vo(Noao  + 11d2«)I 


D 


(65) 


Equation  (62)  then  becomes 


4 
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From  Figure  9 it  is  obvious  that  the  output  of  the  filter  F(s) 
must  be  e/K^^  when  the  input  is  y.  Therefore  F(s)  becomes 


which  is  the  steady  state  optimal  filter  that  results  in  a 
minimum  value  of  the  variance  of  the  error  stated  by  equation 
(64).  (Note  that  the  derived  equations  (64)  and  (68)  agree  with 
equations  (4-105)  and  (4-104)  of  Ref.  13,  pg  63  for  the  case 
where  n=l,  Q=l,  k=l,  Kt=l,  and  C(k)=4R=4A2/aNQ  .) 

p 

Equation  (64)  is  the  miriimum  value  for  due  to  the  effects 
of  fluctuations  of  both  the  process  0 and  the  noise  Nc(t,$). 

The  bandwidth  ag  of  equation  (64)  is  smaller  than  of  equation 
(45)  because  of  the  a pn^ion-i  statistics  of  the  process  0 which 
were  designed  into  the  filter  F(s). 


Suppose  that  a plane  wave  is  incident  on  the  entrance 
aperture  of  a shearing  interferometer.  Assume  that  the 
incident  phasefront  has  some  x-tilt  as  shown  in  Figure  11. 

It  will  be  shown  that  by  varying  the  speed  of  the  rotating 
grating  using  a feedback  network,  the  off  axis  tilt  in  the 
x-direction  can  be  determined.  That  is,  optical  rather  than 
electrical  feedback  is  used. 

The  optical  field  can  be  written  as  (Ref.  18,  pg  49) 

E (x , y ) = EQexp [ jkcosex]  (69) 

From  the  earlier  analysis  of  the  shearing  interferometer  the 
linear  velocity  of  the  grating  was  v = rwg  (a  constant  speed) 
and  the  grating  transmittance  function  g(x^-vt)  was  stated  by 
equation  (3).  If  the  linear  velocity  is  allowed  to  vary,  then 
the  grating  transmittance  function  becomes  g(x^-r/^w^(a)da) 
where  vt  has  been  replaced  by  r/^w(a)da.  The  spatial  Fourier 
transform  of  this  new  grating  function  becomes 

F[g(x1-r/^w1(o)da)  ] « G ( x2~r /^w.^  ( a ) da ) 


Figure  11.  Plane  wave  with  x-tilt 
incident  on  a shearing  interferometer 


kx 

= £ exp[-j^(r/QW1(a)da)  ]An6(x2- ^f)5(fy)  (70) 

The  electric  field  at  the  detector  plane  becomes  (from  equation 


ud(x2,y2)  = w(-x2,-y2)  0G(x2-r/gW1(a)da) 


-l-oo  1/y  ^ f" 

/_oo  F,oexp[jkcos0(x-x'  ) ] exp  [ -j— jr-r/gW1  ( a ) da  ; 


• l A 6(x'-  ^)dx'6lfj 
n n u y 


4 


-j-  CO  1/-  Y t 4" 

+ E0exp[,ikcosQxJ/_co  exp  [ -j— (rf^  ( a )da  - cosO)] 


' h Ani(x'- 


n X "P 

= EQexp  [jkcosGx  J E Anexp  [ jkq^-— ] 6 ( f ) 

y 

where  q = r/gW^(a)da  - cosG  . The  intensity  becomes 


IH  = uau*a  = lEnl<m£  A2  + E A A* 
d d d 1 0 1 m=n  n , n m 


exp  [-jkq-^j(n-m)  ] } 


where  again  n can  only  take  on  the  values  of  0,  +1,  and  -1 
so  that  double  and  higher  order  frequency  terms  can  be 
neglected.  This  results  in  the  following  bandpass  portion  of 
equation  (72) 


!d  = |Eq|2[(I  A2)  + A0A*_1exp[-jk(^)q] 


+ A0A*+1exp[+jk(^p)q] 


+ A_1A*0exp[+jk(^i-)q] 


+ A+1A*0exp[-jk(^)q] 
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where  as  shown  in  Table  1 (Chapter  2)  AQ=.5,  A^  = ( i/ir  )exp  [+j  90 0 ] 
and  A_1=(l/ir)exp[-j90£’ ] . Therefore  equation  (73)  becomes 

Id  = |EQ  | 2{  ( E A2)  + —■[  j cos  ( kq-^p)  + sin(kq^) 

- jcos(kq-^p)  + sin(kq-^p)  - jcos(kq-^p) 

+ sin(kq-p)  + jcos(kq-^p)  + sin(kq-^p-)  ] } 

= |E0|2{(£  An}  + fsin[^r/^w1(a)da  - ^pcos©]  (74) 


2 7T  f*  ^ ^ t 

Let  ip~  d cos G and  Wgt  + ip  = ^pr/ gW^( a )da  then  ignoring  DC 
terms  and  amplitudes,  equation  (74)  becomes  of  the  form 

Id  = sin(wQt  + J-^)  (75) 

This  signal  could  be  used  as  the  input  to  the  feedback  system 
of  Figure  12  where  i ji-\p  represents  an  error  signal  input  to  a 
Kalman  filter  as  described  in  the  last  section.  When  >Jj  = i|> 
then  no  command  will  be  sent  to  the  grating  and  it  will  rotate 
at  a quiescent  velocity,  and  the  off-axis  tilt  will  be  known. 

Therefore  by  employing  optical  feedback  to  the  grating 
and  alternate  technique  for  determining  off-axis  tilt  is 
possible . 
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Coherent  Addition  of  Two  Detector  Arrays 


One  of  the  major  drawbacks  to  using  a shearing  interfer- 
ometer for  phasefront  sensing  or  measurement  is  the  poor 
utilization  of  input  energy.  Much  of  the  input  energy  is 
lost  due  to  the  less  than  unity  transmittance  of  the 
optics  (which  includes  the  lenses  and  Ronchi  grating).  In 

this  section  it  is  shown  that  the  addition  of  another  detector 
array  appropriately  located  can  increase  the  signal  to  noise 
ratio . 

Assume  that  a detector  output  of  a shearing  interferometer 
is  expressed  as 


s^Ct)  = Dsin(wQt  + e(t))  + n-j^Ct)  (76) 

where  n.^  is  a zero  mean  random  process.  The  signal  to  noise 
ratio  for  equation  (76)  is 


SNR1 


D 


2NqW 


(77) 


p 

where  D /2  is  the  signal  power,  NQW  is  the  noise  power,  and  V.’ 
is  the  bandwidth. 

The  arrangement  of  Figure  13  shows  a possible  configuration 
for  the  addition  of  another  detector  array.  The  rotating  Ronchi 
grating  has  been  tilted  at  a slight  angle  so  that  the  reflected 
portion  of  the  field  is  recieved  by  the  additional  detector 
array.  A signal  at  the  second  detector  array  becomes 


s2(t-td)  = Dsin(wQ(t-td)  + e(t-td))  + n2(t)  (78) 
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Figure  13-  Arrangement  for  additional  detector  arra 


where  td  is  the  delay  time  with  respect  to  the  original  detector 
array,  and  n^Ct),  ^(t)  are  uncorrelated  and  have  the  same 
power  of  NqW.  The  delay  time  can  be  represented  as 


t,  = 


width  of  spoke  where  image  is  located 


linear  velocity  of  the  grating 


If  the  output  of  the  original  detector  array  can  be  delayed  by 
td  as  shown  in  Figure  13,  then  the  two  signals  s^(t)  and  S2(t) 
can  be  coherently  added: 
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st(t-td)  = + s2U-td) 


= 2Dsin(wQ(t-td)  + 0(t-td))  + n1(t)  + n^t)  (80) 


For  equation  (80)  the  signal  power  is  2D^  and  the  noise  power 

2 

becomes  E[(nd(t)  + ] = 2NQW.  The  signal  to  noise  ratio 

now  becomes 


snr2  * sb 


(81) 


which 

This 


into 

l 


is  an  improvement  of  3dB  over  the  result  of  equation  (77). 
result  assumes  that  the  entire  input  field  is  reflected 
the  second  detector  array  rather  than  scattered. 


' 


1 
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IV  Representation  of  a Distorted  Phasefront 


A geometric  representation  of  a distorted  phasefront  will 
be  provided  in  this  chapter  (Refs.  19,  PS  1427-1435;  20, 
pg  144-154;  21,  pg  464-466).  This  representation  can  be  made 
up  of  parameters  such  as  tilt,  defocus,  and  astigmatism  which 
are  understood  optics  terms. 

Fried  (Ref.  19)  derives  a description  of  the  phasefront 
shape  by  fitting  the  distorted  phase  profile  with  a polynomial. 
It  is  assumed  that  the  distorted  phase  profile  is  received  by 
a circular  aperture  of  diameter  d.  A set  of  orthonormal 
polynomials  which  are  closely  related  to  Zernike  polynomials 
(Ref.  21)  are  used.  The  six  orthonormal  polynomials  are 
defined  as  follows  (where  R = d/2  and  p defines  the  polar 
representation  of  x and  y coordinates  in  the  receiving  aperture) 


Zx(p)  = (1/ttR2)1/2  (82) 

which  is  used  to  represent  an  overall  phase  shift, 

Z2(p)  = (4/urV/2x  (83) 

and 

Z3(p)  = ( 4/ uR^ )1/2y  (84) 

are  used  to  represent  x and  y tilt  respectively.  The  focus 
term  is  represented  by 

z4(p)  = (12/rrR6)1/2(x2  + y2  - R2/2)  (85) 

and  astigmatism  is  represented  by 

Z5(p)  = ( 6/ uR6 ) 1/2 ( x2  + y2)  (86) 
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and 


Z6(p)  = (24/TrR6)1/2xy  (87) 

The  distorted  phase  can  be  approximated  by 

*(p>  = J,  a Z (p)  (88) 

J J J 

where  the  expansion  coefficients  a^.  are  defined  by 

a,  = /+°°  dp  £ ( n ; d ) (j>  ( p ) Z . ( p ) (89) 

where  n=  |p|  , C(  n ;d)  is  the  receiving  aperture  function 

1 if  | p | <R 

£(n;d)  = (90) 

.0  if  |p|>R 

and 

a^  « degree  of  phase  shift  of  the  phasefront 
a^  “ extent  of  random  x-tilt 
a^  11  extent  of  random  y-tilt 
aj,  <*  extent  of  defocus 

a5  + a6  “ extent  astigmatic  deformation 

As  mentioned  earlier,  the  measurements  made  by  the  shearing 
interferometer  consist  of  twelve  phase  difference  measurements 
0 « 4>(x  -s)  - ^(x^+s)  . The  basic  definition  of  a derivative 
implies  that 


(91) 


9 ~ ix  ^or  fy  for  y_shear) 


when  the  shear  distance  is  small.  Equation  (88)  now  becomes 


(for  x-shear) 


9<t>(x,y) 


6 3Z  lx,y) 

= , Z,  a. — i 

J=1  J 9x 


In  matrix  form,  equation  (92)  becomes 


! = ]a  (93) 

where  0_  is  a (12x1)  vector  of  phase  difference  measurements 
provided  by  the  phasefront  sensor  (six  from  the  x-channel 
interferometer  and  six  from  the  y-channel  interferometer).  The 
quantity  [-rr(Z)]  is  a (12x6)  matrix  of  twelve  rows  of  appropriate 

a X — 

partial  derivatives  of  equations  (82)  to  (87),  and  a is  a (6x1) 
vector  containing  (a^-a^).  Equation  (93)  has  twelve  equations 
and  six  unknowns  (a^-a^).  This  overdetermined  system  can  be 
solved  by  the  method  of  least  squares.  (See  Ref.  22  for  a 
contrasting  approach.)  Since  [— ( Z ) ] is  rectangular,  it  does 
not  posess  an  inverse;  therefore,  in  order  to  find  a the 
technique  of  generalized  inverses  must  be  used.  (A  similar 
procedure  for  generalized  inverses  is  found  in  Refs.  23  and  24.) 
Equation  (93)  can  be  rewritten  as 


a = t^(Z)]  + e 


where 


39 


(95; 


0 

which  is  the  generalized  inverse  of  [-r—  (Z)],  and  a is  a 

oX  — — 

best  approximate  solution  of  through  . 

In  practice  0 is  a (12x1)  vector  of  e's  which  are  obtained 

from  the  phaselocked  loop  of  the  previous  chapter.  Therefore, 

a is  actually  found  from  equation  (9*0.  The  relationship 
2 2 

between  o ; and  a - is  derived  as  follows 
— a 

cov [a ] = E [ ( a - E[a])(a  - E[a])T]  (96) 

T 

= [f^(Z)  ]+E[(e  - e[£]  (§  - e[9])t)  [— (z;  ]+  (97; 

= t|I(Z)]+cov[0][|7(Z) I+T  (98) 

Based  on  earlier  results  the  diagonal  terms  can  be  found.  The 
cross  terms  of  cov[9_]  are  nonzero  since  the  0's  are  correlated. 
These  cross  terms  cannot  be  easily  found. 

In  summary,  the  output  of  the  phasefront  sensor  consists  of 
twelve  phase  difference  measurements  which  can  be  used  to  fin  : 

a in  a computer  implementation.  Then  the  obtained  values  of 
a^  through  can  be  translated  into  control  commands  to  corr-  c* 
for  random  tilting  of  the  phasefront,  defocus,  and  astigmatism 
introduced  by  the  turbulent  atmosphere. 
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V Conclusion 


Summary 

This  paper  first  presented  a description  of  operation  for 
a shearing  interferometer.  It  was  shown  that  the  interferometer ' r 
output  were  signals  with  phases  which  were  a function  of  a 
phasefront's  spatial  phase. 

The  output  signal  processing  was  studied  using  a basic 
demodulation  scheme  and  then  compared  to  a design  using  phase- 
locked  loops  and  some  optimal  linear  filtering  theory.  Some  •] 

basic  phaselocked  loop  theory  was  provided. 

A possible  configuration  using  optical  feedback  for 
determining  off-axis  tilt  of  a plane  was  discussed.  It  was 
found  that  by  varying  the  speed  of  the  rotating  grating,  the 
off-axis  tilt  of  an  incident  planar  phasefront  could  be 
determined. 

An  improvement  in  signal-to-noise  ratio  was  achieved  by 
a slight  modification  of  the  shearing  interferometer.  This 
modification  consisted  of  adding  an  additional  detector 
array  to  collect  energy  reflected  off  of  the  rotating  grating 
and  coherently  adding  the  signals  of  both  detector  arrays. 

Finally,  a geometric  representation  of  a distorted  phase- 
front  was  provided.  It  was  shown  that  the  distorted  phasefront 
could  be  represented  by  terms  such  as  tilt,  defocus,  and 
astigmatism. 
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In  this  paper,  noise  contributed  by  dark  current  and 


thermal  effects  was  ignored.  These  effects  on  the  performance 
of  several  signal  processing  schemes  should  be  studied. 

The  fact  that  the  phase  of  each  signal  is  spatially 
correlated  could  be  used  to  design  a better  signal  processing 
arrangement . 

Finally,  the  study  of  alternate  systems  for  phasefront 
sensing  could  be  contrasted  to  the  operation  of  the  shearing 
interferometer. 

The  problem  of  least  square  fitting  a phasefront  estimate 
to  find  an  array  of  phase  difference  measurements  is  discussed 
in  Ref.  25,  pg  370. 
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Appendix 


Intensity  at  the  Shearing  Interferometer  Detector  Plane 
assuming  Sheared  Amplitudes  are  not  Equal 

As  stated  in  Chapter  2 , the  detector  plane  intensity  can  be 
expressed  as 


^d^X2’^2  ^ ^d ^ x2 ’ ^ 2 ^ ^d ^ x 2 } ^ 2 ^ 


" [V0  + V+1  + V-11[V0  + V + v-i>* 

- V*  + W + W + Vo  + v+iv-i 


+ V5  + V?1  + Vh  + V*1  (A-l) 


where  * denotes  the  complex  conjugate  of  the  designated  quantity. 
The  detector  electronics  filters  out  the  DC  terms  (V^V*  > 

> an^  V 1V*1)  and  the  time  varying  terms  of  2Wq  (V+^V*^ 
and  V Therefore  the  information  bearing  portion  of 

equation  (A-l)  becomes 


^d  ^x2  * ^2  ^ 


= fV  V*  + V V*  1 + TV  V*  + V V*  1 
1 +lv0  0 +1J  1 -1  0 0 -1J 

= 2Real{VQV*1)  + 2Real{V_]V§} 

■ 2A0AhAinc<x2’VAinc(x2+s-V 

• cos[wQt  - 4.(x2+s,y2)  + <j>(x2,y2)] 

+ VRohncV-Vhnc'V’V 

. cos [ w q t + 4> (x2— s ,y2 ) - <f>(x2,y2)] 


(A-2) 
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where  Vq  , V+^  , and  V ^ are  defined  in  equation  (10). 

Neglecting  common  factors  equation  (.A-2)  becomes 

Jd  “ Alnc(x2+S>y2)cos[w0t  ~ ^(x2+s,y2)  + <j>(x2,y2)) 

+ Ainc(x2_s,y2)cOstw0t  + ( x2-s , y2 ) - <|>(x2,y2)] 

( A— 3 ) 


which  is  of  the  form 


A’cos[w0t  +<(>']+  A"cos[wQt  + <j>"  ] (A-4) 

Using  phasor  notation  with  an  implied  radian  frequency  of  w^ 
equation  (A-4)  becomes 


[(Ax  + Ax)2  + (Ay  + A")2]1/2exp[jtan~1(-^-y-  ■*  -j|)] 


( A—  5 ) 


where  A'  = A’coscf>'  , A'  = A'sin^'  , A"  = A"cos<f>" 


and 


A"  = A"sin<|>’’.  Including  the  frequency  dependence  of  the  phasor 

«/ 

expression,  equation  (A-5)  can  be  written  as 


KAi  + Af2  + (A’  + AJ) 

Equation  ( A— 6 ) is  a general 
portion  of  the  intensity  at 
For  the  case  discussed 

Ainc(x2,y2)  = Ainc(x2+S,y2) 


A » + a" 

2]1/2cos[wQt  + tan_1(^-T- ] (A-6) 

X X 

result  for  the  information  bearing 
the  detector  plane. 

in  Chapter  2 where  it  is  assumed  that 
“ Ainc(x2_s,y2)  fchen 
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(A- 7) 
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